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The  Pioneering  Research  Laboratory  on  Physics  of  Fine  Particles  was  established  in  1962,  at 
the  Ohio  Agricultural  Experiment  Station,  Wooster,  by  the  Agricultural  Research  Service  of  the 
United  States  Department  of  Agriculture,  in  recognition  of  the  need  for  continued  fundamental  re- 
search on  fine  particle  behavior  and,  in  particular,  the  criticahty  of  such  studies  to  agriculture. 

In  agriculture,  pesticides  in  both  hquid  droplet  and  sohd  particle  form  are  used  extensively,  but 
their  use  is  frequently  hampered  by  lack  of  optimum  efficiency  and  precision  of  application,  despite 
the  best  efforts  of  manufacturers  and  applicators.  This  results  in  an  excessive  economic  burden  on 
agriculture  in  the  maintenance  of  the  best  of  quahty  in  agricultural  produce,  as  well  as  an  increased 
chance  of  contaminating  surrounding  areas.  Problems  of  air  pollution  are  also  of  prime  pertinence 
and  concern.  It  is  desired  that  the  results  of  these  investigations  in  fine  particle  behavior  and  its 
inherent  subject  matter  will  serve  agriculture  and  its  allied  industries,  and  other  industries  with 
similar  vexations,  in  the  alleviation  of  these  problems. 
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PARTICLE-SIZE  DISTRIBUTIONS  IN  BROWNIAN  MOTION 


Ross  D.  Brazee  ^ 


A method  of  incorporating  a particle-size  distribution  into  statistical  treatments  of  Brownian 
motion  theory  is  considered.  Probabilistic  statements  for  the  joint  distribution  of  position  and 
particle-size  and  conditional  distributions  of  position  on  particle-size  and  particle-size  on  posi- 
tion, as  functions  of  time  are  developed.  Chandrasekhar’s  problem  on  the  effect  of  Brownian 
motion  on  gravity  sedimentation  is  used  as  an  illustration  for  determination  of  a joint  distribution. 


INTRODUCTION 


The  apphcations  of  probability  theory  to  the  dynamics  of  Brownian  motion  are  familiar.  In  this  note 
we  will  consider  a method  of  expressing  the  existence  of  a particle-size  distribution  in  such  problems. 
For  a thorough  discussion  of  the  probabilistic  concepts  and  notation  used  in  this  paper,  the  reader  will  find 
Feller  - and  Cramer^  {pp.  267—70,  291-92)  helpful. 


PRELIMINARY  CONSIDERATIONS 


First  imagine  that  we  have  N particles  which  may  be  put  into  n discrete  classes  of  characteristic  size 
Sj{j=\,  2,  n)  with  Nj  particles  in  sj.  The  particles  are  then  subjected  to  an  unspecified  process  in 

time  which  continually  rearranges  them  among  M cells.  Each  of  the  p cells  is  assigned  an  identification 
number  i—l,2,...,  M,  and  we  specify  that  the  total  number  of  particles  shall  remain  constant. 

At  time  t,  we  stop  the  process  and  perform  an  accounting  of  the  particle  arrangement.  We  count  the 
total  number  Nj  of  particles  of  all  sizes  in  cell  z.  The  number  Njj  of  particles  of  size  s,  in  cell  z is  then 
counted.  This  will  give  the  necessary  information  to  estimate  several  probabilities  which  are  required, 
since  we  could  imagine  the  counts  to  be  arithmetic  mean  values  of  several  trials. 

The  probability  of  finding  any  particles  in  cell  i,  at  time  t,  irrespective  of  their  size  sj  will  be 


= (1) 

The  probability  of  finding,  at  time  t,  any  particles  of  size  sj  in  cell  i will  be 

(2) 

The  probability  of  finding  any  particles  of  size  sj  in  the  entire  system  will  be  independent  of  time  and 
given  by  the  particle-size  probability  function 


(3) 


* Pioneering  Research  Laboratory  on  Physics  of  Fine  Particles,  Agricultural  Engineering  Research  Division,  ARS,  USDA,  Wooster, 

^ W.  Feller.  An  Introduction  to  Probability  Theory  and  Its  Applications.  2d  ed.,  Vol.  i.  Chap.  v.  New  York.  1957. 

3 H.  Cramer.  Mathematical  Methods  of  Statistics.  575  pp.  Princeton,  N.J.  1946. 
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The  conditional  probability  that  a particle  will  be  found  in  cell  i at  time  t on  the  hypothesis  that  it  is  of 
size  Sj  will  be  given  by 


p{i\j-  111 . 


(4) 


The  probability  that  a particle  of  size  Sj  will  be  found  at  time  t on  the  condition  that  we  search  in  cell  i 
will  be 


(5) 


Ni 

Eq.  (5)  is  a statement  of  the  particle-size  distribution  in  space  as  a function  of  time. 

In  (1),  (2),  and  (4),  and  (5)  the  variable  t is  not  distributed,  but  the  probabilities  are  functions  of  time. 


THE  DISCRETE  CASE  IN  ORDINARY  SPACE 


The  development  is  similar  to  that  of  the  preceding  section  if  the  cells  are  replaced  by  M discrete 
regions  in  ordinary  space.  In  this  case,  the  location  of  each  region  is  specified  in  terms  of  a position 
vector  r;  with  components  of  scalar  magnitude  rn,  rj2,  and  ns. 

The  conditional  probabihty,  at  time  t,  of  the  event  ^ = rj  relative  to  the  hypothesis  7]  = Sj  will  be  as  in 
(4),  viz. 


P{^  = Ti\r]^Sj;t) 


P{^  = T^u  V = Sj;  t) 
P{r)  = Sj) 


^Pij 
P.j  ’ 


(6) 


where  pij  is  the  joint  probabihty  of  the  events  i = Ti  and  r]  = Sj,  and  p.j  is  the  marginal  probability 


IVI 

Pj~^  Pij- 


(7) 


The  conditional  probabihty,  at  time  t,  of  the  event  r]  = Sj  on  the  hypothesis  that  ^ — wiU  be  given  by 


P{v  = Sj\^  = ri;t) 


P{^  = T^i,V  = Sj;t) 
P{^  = ri;t) 


(8) 


as  in  (5).  Since  the  events  r)  = sj  are  mutually  exclusive,  we  may  write 


P{^^ri;t)  = '^[P{^  = ri\p  = Sj:t)P{7)  = Sj)] 
j 

n n 

~ Pi-  ~ ^ Pij  ~ ^ Pi\jP  j- 
j j 


(9) 


Hence,  (8)  becomes 


P{ri  = Sj\i  = ri;t) 


P{i  = ri,7]  = sj;  t) 

^ [P(^  = 'ri\r}  = Sj;  t)P{r}  = Sj)] 


Pij  Pa . 

71  . 

V P^- 

^Pi\jP-j 


(10) 
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THE  CONTINUOUS  CASE  IN  ORDINARY  SPACE 


Let  f{r,  s;  t)  be  the  joint  frequency  function  in  time  for  the  variables  r and  5.  Then  conditional  proba- 
bilities may  be  written  for  the  continuous  case  in  a manner  analogous  to  the  discrete  systems.  We  assume 
that  the  distribution  function  F(r,  which  corresponds  to  the  probabihty  function  P{^  ^ s;  t), 

viz, 


F{ri,  T2,  rs,s;t)  =P(^i  ^ 

ri,  ^2  ^ T2,  ^3  ^ r3,  T7  ^ 5,- 1)  , 

(11) 

is  everywhere  continuous. 

and  that 

C/5 

II 

_a^F(n,  t2,  t3,  s;  t) 
dridr2drsds 

(12) 

is  continuous  except  for  a 

finite  number  of  curves. 

Then  for  any  set  S,  it  follows  that 

PiS,  t)  = [ 

‘'S 

y(r,  5,-  t)dvds, 

(13) 

and  in  particular. 

[”/(r:,r2 

J 0 J— 00  J —CO  J —CO 

, rs,  s;  t) dr\dr2dr^ds  = 1. 

(14) 

The  marginal  distribution 

of  ^ will  have  the  probability  function 

P{^^ 

r;t)=h  h h [7(u, 

J —00  J —00  J —00  J 0 

v;t)dvd\i=  f^^/i(u;  t)du, 

J —00  J —00  J —y: 

(15) 

where 

/i(r;0=j 

f{r,s;t)ds. 

0 

(16) 

Similarly,  the  marginal  distribution  of  r\  will  possess  the  probability  function 

noo  r 00  r 00  c s 

f{\x,v;t)d\xdv=  \ f2{v)dv, 

-X*/— ooj— 00  Jo 

where 

r X r X r X 

f2(s)=  f(r,s;t)dr 

J —00  j —00  •'—X 

is  precisely  the  overall  particle-size  distribution  frequency  function. 

The  conditional  probability  of  the  event  ^ r relative  to  the  hypothesis  s < iq  < s~\-  As  is 


(17) 


(18) 


P{^  ^r|5<T7<5  + As;  t) 


_ .P(^  r,  5 < 17  < 5 + A5;  0 _ J 

Cs+As 

[ J 

L7..J 

C^i 

1 f{Y,s;t)drds 

— X 

P{s<  iq<  s-\-  As)  j 

I"  s+As  ^ 

— X J —X  J 

1 /(r,  s;  t)drds 

— 00 

If  As  tends  to  zero,  then 


(19) 
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(20) 


lim 

As 


0 


P(^  ^ r|5  <7]  < s -\-As;  t)— 


p r r fir,  s:t)di 

J —cc  — ac  d — oe 

rx  roo  roc 

I f(r,s;t)di 

J —rn  •'—as  —cc 


p p rf(r,s;t)dr 

r —CC  J —PC  J —PC 

f2(s) 


Eq.  (20)  is  the  conditional  distribution  function  in  r for  fixed  s.  If /(r,  s;  t)  is  continuous  in  r,  it  may  be 
differentiated  with  respect  to  the  components  of  r to  obtain  the  conditional  frequency  function 

/(rk0=^^-  (21) 

The  probability  of  the  event  rj  ^ s on  the  hypothesis  r < ^ < r + Ar  is 


r„+Ar„  ,.r„+Ar„  -r, +Ar, 


P(i7^5|r<^<  r + Ar;  t) 


P{r  <i  < r + Ar;  ^) 
If  Ar  is  permitted  to  approach  zero,  (22)  will  yield 


r'3  3 r'2^  2 r 11  r” 

P(T7^5,r<^<r  + Ar;0  _Jr3  j j j^f(r,s;  t)dsdr 

r ro+Aro  r r„+Ar„  f r.+Ar,  f x 

r r f(r,s;t)dsdr 
3 J ^2  J ^1  Jo 


(22) 


lim 

Ar—^0 


P(r)  ^ s\r  < r + Ar;  t)  = 


^ S ^ S 

f(r,s;t)ds  j f{r,s;t)ds 
J n Jq 


f{T,s;t)ds  fiir;  t) 


(23) 


which  is  the  conditional  distribution  function  in  5 for  fixed  r.  If  /(r,  s;  t)  is  continuous  in  5,  it  may  be 
differentiated  with  respect  to  s to  determine  the  conditional  frequency  function 

/(s|r;0=^7^  (24) 

/i(r;  t) 

_ f{r\s;t)f2is) 

[ f(rjs;  t)/2(s)ds 
J 0 

Eq.  (24)  is  an  expression  for  the  frequency  function  of  the  particle-size  distribution  at  point  r at  time  t, 
and  is  analogous  to  (10)  for  the  discrete  case. 


AN  EXAMPLE  OF  DETERMINATION  OF  THE  JOINT  FREQUENCY  FUNCTION 


Chandrasekhar  ^ considers  a one-dimensional  problem  on  the  effect  of  gravity  on  Brownian  motion. 
We  wiU  consider  this  problem  further  in  the  fight  of  the  foregoing  developments  by  indicating  the  deter- 
mination of  the  joint  frequency  function.  We  remark  that  so  far  the  formalism  has  not  been  specifically 
limited  to  Brownian  motion. 

Chandrasekhar  applies  the  Smoluchowski  equation, 

V'(9/3-2Vm;  — K^-^zc),  (25) 


**  S.  Chandrasekhar,  Revs.  Modern  Phys.,  15:  57—9.  1943. 
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where  w is  the  position  frequency  function  for  \ and 


with 


II 

(26) 

^ m 

(27) 

In  (26)  and  (27)  k is  Boltzmann’s  constant,  T,  is  the  absolute  temperature,  ^ is  the  coefficient  of  viscosity 
of  the  surrounding  medium,  and  s and  m are  the  characteristic  radius  and  mass  of  the  particle.  Stokes’ 
law  viscous  dissipation  is  assumed.  The  vector  K represents  the  acceleration  caused  by  an  external 
force  field.  The  semi-infinite  region  z > 0 in  a stiff  medium  is  treated,  with  the  positive  z-axis  in  the 
vertical  direction  and  gravity  as  the  external  field,  so  that 


Kx  = ^',  Ky  = 0\  Kz=^—  [l-{po/p)]g. 


(28) 


where  p and  po(=^  p)  are  the  densities  of  the  particle  and  surrounding  medium,  respectively,  and  g is  the 
acceleration  due  to  gravity.  In  this  situation,  (25)  becomes 


Po\gdiv 


dt  ^ 


\ p / 13  dz  ^ 


(29) 


which  shows  that  the  system  may  be  considered  as  a diffusion  process.  Chandrasekhar  writes 
so  that  (28)  takes  the  form 

rtui  riu) 

(31) 


dw  p.  d^w  . dw 


and  then  imposes  the  initial  and  boundary  conditions 


lim  w = d{z  — zo); 
t^O 

D — -\-cw  = 0,  (z  = 0,  t > 0). 
dz 


(32) 


(33) 


where  8(z  — zo)  is  the  Dirac  delta  function  which  indicates  that  all  the  particles  are  initially  at  z = zq.  The 
solution  of  (31)  subiect  to  the  conditions  (32)  and  (33)  is 


w{z,  t;  Zo)  - 


1 


2{TrDt)'^ 


exp 


■ (z  + zo)“/4Z)^ 


+ exp 


— (z  + zoTI^Dt 


exp 


c 

- ^(2-zo)-^^ 


2D  \2{Dt)'^  ) 


(34) 


Chandrasekhar  makes  several  interesting  observations  about  the  physical  implications  of  the  solution  (34) 
which  we  shall  not  discuss  here.® 

In  (34),  the  constants  c and  D contain  s.  Hence,  in  a system  having  a particle-size  distribution, 
(34)  may  be  written 


w{z,  t;  Zo)  = w(z,  5,  t;  zo).  (35) 

In  applying  (34)  to  a many-particle  system,  it  is  presumed  that  each  particle  obeys  it  independently  of  the 
other  particles.  Similarly,  we  now  assume  that  in  the  size-distributed  system,  the  particles  behave  inde- 

^ See  pp.  58—59  of  footnote  4. 
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pendently,  so  that  (34)  may  be  considered  as  the  conditional  frequency  function.  Consequently,  we  now 
write 


f{z  I5;  ^)  = w{z,  t;  zo). 

By  virtue  of  (21),  the  joint  frequency  function  is 

f{z,  s;  t)  = w{z,  t;  zo)f2{s). 

If,  e.g.,  the  particle-size  is  log-normally  distributed,  i.e.. 


f2{s) 


1 


{27t)'^  s{lno-) 


exp 


{Ins  — ln{s)y 


2{lncry 


(36) 


(37) 


(38) 


where  (s)  is  the  geometric  mean  particle  size  and  cr  is  the  geometric  standard  deviation,  then  combination 
of  (37),  (38),  and  (34)  yields 


/(2,  5;  t) 


1 


27T{2Dt)'^  s{lncr)  [ 


exp  [—  (z  — zo)^/4<Dt]  + exp  [ — (z  + zo)^/4Z)^] 

^z  + zo—  ct 


c{TTt I D)'^  erfc 


2{Dtf 


exp 


exp 


c 


(Ins  — ln{s)y 


2{lnarY 


(39) 


Eq.  (39)  cannot  be  written  in  the  form 


/(z,  5;  t)  = </)i(z,  t)4>2{s,  t),  (40) 

which  points  out  that  z and  s are  not  independent.  At  least  formally,  (39)  could  be  used  to  obtain  an  ex- 
pression for/(5,  z;  t),  the  particle-size  frequency  function  for  point  z at  time  t. 
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